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ABSTRACT 

We consider the influence of microlensing on different spectral bands of lensed QSOs. 
We assumed that the emitting X-ray, UV and optical regions are different in size, but 
that the continuum emission in these spectral bands is originating from an accretion 
disc. Estimations of the time scales for microlensing and flux amplification in different 
bands are given. We found that the microlensing duration should be shorter in the 
X-ray (several months) than in UV/optical emitting region (several years). This result 
indicates that monitoring of the X-ray variations in lensed QSOs that show a 'flux 
anomaly' can clarify the source of this anomaly. 

Key words: accretion, accretion discs - gravitational lensing - galaxies: active - 
ultraviolet: galaxies - X-rays: galaxies 



1 INTRODUCTION 



Recent observational and theoretical studies suggest that 
gravitational microlensing can induce variability not only 
. in opti cal light, but also in the X-ray emission of lensed 
QSOs (IChartas et al.1 I2OO2I 12004 Ibai et al. I I2003I . |2004 
I Jovanovicll2006l ; IPopovic et al . 2001, 2003a b, 2006a b). Vari- 



ability studies of QSOs indicate that the size of the X- 
ray emitting region is significantly smaller (~ several light 
hours) , than the optical and UV emitting regions (~ several 
light days). 

Gravitational lensing is achromatic (the deflection an- 
gle of a light ray does not depend on its wavelength), but it 
is clear that if the geometries of the emitting regions at dif- 
ferent wavelengths are different then chromatic effects could 
occur. For example, if the microle ns is a binary star or if 



the microlensed source is exten d ed (iGriest fc Huf l992. 



1993; 



Bogdanov and Chcrepashchuk Il995allbt [ Zakharov 1 19971 : 
Zakharov fc Sazhinl Il998l ; iPopovic fc Chartasl l2005h differ- 
ent amplifications in different spectral bands can be present. 
Studies aiming to determine the influence of microlensing 
on the spectra of lensed QSOs need to take into account 
the complex structure of the QSO central emitting region 
l|Popovic fc Cha rtas 2005). Since the sizes of the emitting 
regions are wavelength dependent, microlensing by stars in 
the lens galaxy may lead to a wavelength dependent mag- 
nification. For example, Blackburne et al. (2006) reported 
such a 'flux anomaly' in quadruply imaged quasar 1RXS 



J1131-1231. In particular, they found discrepancies between 
the X-ray and optical flux ratio anomalies. Such anomalies 
in the different spectral band flux ratios can be attributed 
to micro- or milli-lensing in the massive lensing halo. In 
the case of milli-lensing one can infer the nature of sub- 
structure in the lensing galaxy, which can be connected to 
Cold Dark Matter (CDM) structures (see e.g. Dobler & Kee- 
ton 2006). Besides microlensing, there are several mecha- 
nisms which can produce flux anomalies, such as extinction 
and intrinsic variab i lity. T hese anomalies were discussed in 
IPopovic fc Chartasl (|2005^ where the authors gave a method 
that can aid in distinguishing between variations produced 
by microlensing from ones resulting from other effects. In 
this paper we discuss consequences of variations due to grav- 
itational microlensing, in which the different geometries and 
dimensions of emitting regions of different spectral bands 
are considered. 

The influence of microlensing on QSO spectra emit- 
ted from their accretion discs in the range from the X-ray 
to the optical spectral band is analyzed. Moreover, assum- 
ing different sizes of the emitting regions, we investigate 
the microlensing time scales for those regions, as well as a 
time dependent response in amplification of different spec- 
tral bands due to microlensing. Also, we give the estimates 
of microlensing time scales for a sample of lensed QSOs. 

In Section 2, we describe our model of the quasar emit- 
ting regions and a model of the micro-lens. In Section 3 we 
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discuss the time scales of microlensing and in Section 4 we 
present our results. Finally, in Section 5, we draw our con- 
clusions. 



2 A MODEL OF QSO EMITTING REGIONS 
AND MICROLENS 

2.1 A model of the QSO emitting regions 

In our models we adopt a disc geometry for the emit- 
ting regions of Active Galactic Nuclei (AGN) since the 
most widely accepted paradigm for AGN incl udes a super- 
massiv e black hole fed by an accretion disc. iFabian et al.l 
(1989) calculated spectral line profiles for radiation emitted 
from the inner parts of accretion discs and later on such 
featur es of Fe Ka lines were discovered by iTanaka et al.l 
l| 19951 ) in Japanese ASCA satellite data for Seyfert galaxy 
MGC-6-30-15. Moreover, the assumption of a disc geome- 
try for the distribution of the emitters in the central part 
is support ed by the spe c tral s hape of the Fe Ka line in 



tions dZakharov & Remnl 199S 


, 20021; Zakharov et all 20031; 


Zakharov & ReDinl 12003 Jblld. 


2004. 120051. 120061: IZakharovl 


2004!. 120071")). On the other hand, verv often a buniD in the 



UV band is present in the spectra of AGN, that indicates 
that the UV and optical continuum originates in an accre- 
tion disk. 

We should note here that probably most of the X-ray 
emission in the 1-10 keV energy range originates from in- 
verse Compton scattering of photons from the disc by elec- 
trons in a tenuous hot corona. Proposed geometries of the 
hot corona of AGN include a spherical corona sandwiching 
the disc and a patchy corona made of a few com pact re- 
ions covering a small fraction of the disc (see e.g. iMalzad 
20071 )1. On the other hand, it is known that part of the 
accretion disc that emits in the 1-10 keV rest-frame band 
(e.g. the region that emits the continuum Compton reflec- 
tion component and the fluorescent emission lines) is very 
compact and may contribute to X-ray variability in this en- 
ergy range. In order to study the microlensing time scales, 
one should consider the dimensions of the X-ray emitting 
region which are very important for the integral flux vari- 
ations due to microlensing. The geometry of the emitting 
regions adopted in microlensing models will affect the sim- 
ulated spectra of the continuum and spec tral line emission 
l|Popovic fe Chartasll2005l ; |jovanovic]|2006l ) . In this paper we 
assume that the AGN emission from the optical to the X-ray 
band originates from different parts of the accretion disc. 

Also we assume that we have a stratification in the 
disc, where the innermost part emits the X-ray radiation 
and outer parts the UV and optical continuum emission. 
To study the effects of microlensing on a compact accre- 
tion disc we use the ray-tracing method, considering only 
those photon trajectories that reach the sky plane at a given 
observer angle (see, e.g. Popovic et al. 2003a,b and refer- 
ences therein). The amplified brightness with amplification 
A(X, Y) for the continuum is given by 



Ic(X, Y\ Eobs) — Ip{E bs,T(X, Y)) ■ A(X, Y), (1) 

where T(X, Y) is the temperature, X and Y axe the impact 
parameters which describe the apparent position of each 



point of the accretion disc image on the celestial sphere as 
seen by an observer at infinity, E b s is the observed energy, 

Ip is an emissivity function. 

In the standard Shakura-Sunyaev disc model (|l973T ). ac- 
cretion occurs via an optically thick and geometrically thin 
disc. The effective optical depth in the disc is very high and 
photons are close to the thermal equilibrium with electrons. 
The surface temperature is a function of disc parameters and 
results in the multicolor black body spectrum. This compo- 
nent is thought to explain the 'blue bump' in AGN and the 
soft X-ray emission in galactic black holes. Although the 
standard Shakura-Sunyaev disc model does not predict the 
power-law X-ray emission observed in all sub-Eddington ac- 
creting black holes, the power law f or the X-ray emissivi ty 
in AGN is usually adopted (see e.g. iNandra et al.l |l999)). 
But one can not exclude other approximations for emissivity 
laws, such as black-body or modified black-body emissivity 
laws. Moreover, we will assume that the emissivity law is 
the same through the whole disc. Therefore we used the 
blac k-body radiation l aw. The disc emissivity is given as 
(e.g. lJaroszvnskil (|l992l )'): 



Ip(X,Y;E) = B[E,T s (X,Y)], 



where 



B(E,T S (X,Y)) 



2E' A 



1 



/j2 c 2 e E/kT s (X,Y) 



(2) 



where c is the speed of light, h is the Planck constant, 
k is the Boltzmann constant and T s (X, Y) is the sur- 
face temperature. Concern ing the standard accretion disc 
llShakura fe Sunavcv] Il973h . here we assumed that (see 
iPopovic et all (|2006al )l: 

T S (X, Y) ~ r- 3/2 {X, Y)(l - r~ 1/2 {X, Y)) 4/5 K, (3) 

taking that an effective temperature in the innermost part 
(where X-ray is originated) is in an interval from 10 7 K to 
10 s K. 



2.2 A model for microlensing 

To explain the observed microlensing events in quasars, one 
can use different microlensing models. The simplest ap- 
proximation is a point-like microlens, where microlensing 
is caused by some compact isolated object (e.g. by a star). 
Such a microlens is characterized by its Einstein Ring Radius 

in the lens plane: ERR = \j < ~'I n ^ >l ^ la or by the corre- 



al 



sponding projection in the source plane: Re — -j=j-ERR = 

— t-, la , where G is the gravitational constant, c is 
c 2 D t 

the speed of light, m is the microlens mass and Di, D s 
and Di s are the cosmological angular distances between 
observer-lens, observer-source and lens-source, respectively. 
In most cases we can not simply consider that microlensing 
is caused by an isolated compact object but we must take 
into account that the micro-deflector is located in an ex- 
tended object (typically, the lens galaxy). Therefore, when 
the size of the Einstein ring radius projection Re of the 
microlens is larger than the size of the accretion disc and 
when a number of microlenses form a caustic net, one can 
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Figure 1. Magnification map of the QSO 2237+0305A image (left) and of a "typical" lens system (right). The white solid lines represent 
three analyzed paths of an accretion disc center: horizontal (y = 0), diagonal (y = —x) and vertical (x = 0). 



describe the microlensing in terms of the crossing of the disc 
by a straight fold caustic (Schneider et al. 1992). The am- 
plification at a point of an exte nded source (accretion d isc) 
close to the caustic is given bv lChangfeRefsdall dHH) (a 
more general expressio n for a magnification near a cusp type 
singul arity is given by ISchneider fc Weiss! (|l992T ) ; IZakharovl 
(199!)): 



A[X, Y) = Ao + K 



^caustic 



(4) 



where Ao is the amplification outside the caustic and K = 
Aof3 is the caustic amplification factor, where /3 is constant 
of order of unity (e.g. Witt et al. 1993). The "caustic size" 
)"caustic is the distance in a direction perpendicular to the 
caustic for which the caustic amplification is 1, and therefore 
this parameter defines a typical linear scale for the caustic in 
the direction perpendicular to the caustic. £ is the distance 
perpendicular to the caustic in gravitational radii units and 
£ c is the minimum distance from the disc center to the caus- 
tic. Thus, 



tea = 



Y_ 



and 



C = Cc + 



(X - Xejtg^ + Y c - Y 



(5) 



(6) 



(7) 



where <f) = a + tt/2. JJ(k(£ — f c )) is the Heaviside function, 
//(k(£ — f c )) = 1, for — £ c ) > 0, otherwise it is 0. k 
is ±1, it depends on the direction of caustic motion; if the 
direction of the caustic motion is from approaching side of 
the disc k = — 1, otherwise it is +1. Also, in the special case 



of caustic crossing perpendicular to the rotating axis k = +1 
for direction of caustic motion from — Y to +Y, otherwise it 
is — 1. A microlensing event where a caustic crosses over an 
emission region can be described in the following way: before 
the caustic reaches the emission region, the amplification is 
equal to Ao because the Heavisied function of equation (4) is 
zero. Just as the caustic begins to cross the emitting region 
the amplification rises rapidly and then decays gradually 
towards Ao as the source moves away from the caustic-fold. 

Moreover, for the specific event one can model the caus- 
tic shape to obtain different parameters (see e.g. Abajas 
et al. 2005, Kochanek 2004 for the case of Q2237+0305). 
In order to apply an appropriate microlens model, ad- 
ditionally we will consider a standard microlensing mag- 
nification pattern for the Q2237+0305A image (Fig. 1 
left). Fo r generating this map we used the ray-shootin 
method dKavser et al.lll986l : ISchneider fc Weisslll98A Il98 
IWambsganss et al.lll990al lbl). In this method the input pa- 
rameters are the average surface mass density a, shear 7 
and width of the microlensing magnification map expressed 
in units of the Einstein ring radius (defined for one solar 
mass in the lens plane). 

First, we generate a random star field in the lens plane 
with use of the parameter a. After that, we solve the Poisson 
equation V 2 ip = 2a in the lens plane numerically, so we can 
determine the lens potential ip in every point of the grid in 
the lens plane. To solve the Poisson equation numerically 
one has to write its finite difference form: 



ipi+i,i + i>i-i,j + ipi,j+l + 4>i,j-i - 4tpij = 2CTi,j. (8) 

Here we used the standard 5-point formula for the two- 
dimensional Laplacian. Next step is inversion of the equation 
((8| using Fourier transform. After some transformations we 
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obtain: 



2 (cos 



2)' 



(9) 



where Ni and N2 are dimensions of the grid in the lens 
plane. Now, using the finite difference technique, we can 
compute the deflection angle a = X7ip in each point of the 
grid in the lens plane. After computing deflection angle, we 
can map the regular grid of points in the lens plane, via lens 
equation, onto the source plane. These light rays are then 
collected in pixels in the source plane, and the number of 
rays in one pixel is proportional to the magnification due to 
microlensing at this point in the source plane. 

Typically, for calculations of microlensing in gravita- 
tionally lensed systems one can consider cases where di- 
mensionless surface m ass density a is some fractio n of f 
e.g. 0.2, 0.4, 0.6, 0.8 l|Trever fc Wambsganssl |2004| ), cases 
without external shear (7 — 0) and cases with 7 = a for 
isothermal sphere model for the lensing galaxy. In this ar- 
ticle we assume the values of these parameters within the 
range generally adopted by o ther authors, in particular by 
ITrever fc Wambsganssl l|2004 ). 

Microlensing magnification pattern for the 
Q2237+0305A image (Fig. I left) with 16 R E on a 
side (where Re ~ 5867 R g ) is calculated usin g the following 
parame ters: a = 0.36 and 7 = 0.4 (see IPopovic et al.l 
(|2006al ). Fig. 2), the mass of microlens is taken to be O.3M 
and we assume a flat cosmological model with SI = 0.3 and 
H = 75 km s _1 Mpc -1 . 

We also calculated microlensing magnification pattern 
for a "typical" lens system (Fig. I right), where the redshifts 
of microlens and source are: zi — 0.5 and z s = 2. In this case, 
the microlens parameters are taken arbitrary: a = 0.45 and 
7 = 0.3 and the size of obtained microlensing pattern is also 
16 R E x 16 R E , where R E fa 3107 R„. 



3 TYPICAL TIME SCALES FOR 
MICROLENSING 

Typical scales for microlensing are discussed not only in 
book s on gravita ti onal l en sing llSchneider. Ehlers fc Falco 
19921 : IZakharovl Il997l ; iPetters. Levine fc Wambsganss 
but in recent pape rs also (see. 



2001) , but in recent pape rs also (see, for example, 
Trever fc Wambsganssl |2004)). In this paper we dis- 
cuss microlenses located in gravitational macrolenses 
(stars in lensing galaxies), si nce optical depth for m i- 



crolensing is then the highest ( Wvithe fc TurneH l2002al lbl: 



I Zakharov. Popovic fc Jovanovic! |2004 120051 ) m comparison 
with other possible locations of gravitational microlenses, 
as for example stars sit uated in galactic cluster s and 
extragalactic dark halos llTadros. Warren fc Hewettl 1 19981 : 
lTotanill2003, ; llnoue fc Chiball2003j ). 

Assuming the concordance cosmological model with 
fltot = 1, n ma ttcr =0.3 and f^A = 0.7 we recall that typ- 
ical l ength scale for microlensing is l|Trever fc Wambsganssl 
l2004h : 



Re = ^2r s ■ ££± « 3.2 • 10 16 ^^5°' 5 cm, (10) 

where "typical" microlens and source s redshifts are assumed 
to be zi — 0.5, z.s — 2. (similar to ITrever fc Wambsganssl 



»), r a = 

c A 

ing to microlens mass m, hrs - 
dimensionless Hubble constant. 

The corresponding 

l|Trever fc Wambsganss| [2004) 



is the Schwarzschild radius correspond- 
H /((75 km/sec)/Mpc) is 

angular scale is 



Re 



2.2 • 10" 



— h-° ! 



(11) 



Using the leng th scale 1|10|) and vel o city scale (say 
v± ~ 600 km/sec as ITrever fc W ambsganss (2004) did), one 
could calculate the standard time scale corresponding to the 
scale to cross the projected Einstein radius 



t E 



(1 + zi 



Re 



25 



m 

w 



-1 7-0.5 
^600 "75 



years, 



(12) 



where a relative transverse velocity «6oo = «x/(600 km/sec). 
The time scale ts, corresponding to a point- mass lens and to 
a small source (compared to the projected Einstein radius of 
the lens), could be used if microlenses are distributed freely 
at cosmological distances and if each Einstein angle is lo- 
cated far enough from another one. However, the estimation 
(|12|) gives long and most likely overestimated time scales 
especially for gravitationally lensed systems. Thus we must 
apply another microlens model to estimate time scales. 

For a simple caustic model, such as one that considers 
a straight fold causticQ, there are two time scales depending 
either on the "caustic size" (r caus ti c ) or the source radius 
(^source). In the case when source radius is larger or at least 
close to the "caustic size" {R B ource > ^caustic), the relevant 



time scale is the "crossing time" I Trever fc Wambsganssl 
l2004h : 



v±{D a /Di 



0.69 Ris w 60 o yj^J h 7°' 5 y ears 
251 Ri 5 Vqo /175 ' 5 days, 



(13) 



where Di and D a correspond to zi — 0.5 and z s = 2, re- 
spectively and R15 = Tisource/10 1 cm. As a matter of fact, 
the velocity perpendicular to the straight fold caustic char- 
acterizes the time scale and it is equal to v± sin (3 where 
is the angle between the caustic and the velocity v± in 
the lens plane, but in our rough estimates we can omit fac- 
tor sin (3 which is about unity. However, if the source ra- 
dius -Rsourcc is much smaller than the "caustic size" r ca ustic 
(^source <S ^caustic), one could use the "caustic time", i.e. 
the time when the source is located in the area near the 
caustic: 



tc 



v±{D a /Di 



D a 



0.69 ris «600 h ^' 5 y ears 

251 ri 5 Vg 00 /if 5 015 days, 



(14) 



1 We use the following approximation for the extra magnification 
near the caustic: 11 - 



£ > £ c where J is the perpen- 



dicular direction to the caustic fold (it is obtained from Eq. (J4j 
assuming that factor K is about unity). 
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Figure 2. The variations of normalized total continuum flux in 
optical (3500-7000 A), UV (1000-3500 A) and X (1.24-12.4 A 
i.e. 1-10 KeV) band due to microlensing by a caustic crossing 
along y = —x direction in the case of Schwarzschild metric. Time 
scale corresponds to "typical 1 ' redshifts of microlens and source: 
zi = 0.5 and z s = 2. The parameters of the caustic are: Aq=1, 
/3=1j k = and its "size" is 9000 R g . Negative distances and 
times correspond to the approaching side, and positive to the 
receding side of accretion disc. In this case, due to k = +1, caustic 
motion is from the receding towards the approaching side (i.e. 
from the right to the left). The source mass is 10 s Mq. The radii 
of optical emitting region are: Ri n = 100 R g , Rout — 2000 R g , 
for UV emitting region: Ri n = 100 R g , R ou t = 1000 R g and for 
X emitting region: R in = R ms , Rout = 80 R g . 
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Figure 3. The same as in Fig. 2. but for r" cauat i c = 2000 R g . 

where ns = r C austic/10 15 cm. 

Therefore, t CVOBB could be used as a lower limit for typical 
time scales in the case of a simple caustic microlens model. 
From equations (|13p and (|14p it is clear that one cannot 
unambiguously infer the source size -Rsourcc from variability 
measurements alone, without making some further assump- 
tions. In general, however, we expect that t cross corresponds 
to smaller amplitude variations than Caustic, since in the 
first case only a fraction of a source is significantly ampli- 
fied by a caustic (due to assumption that -R SOU rcc > ^caustic), 
while in the second case it is likely that the entire source 
could be strongly affected by caustic amplification (due to 
assumption that -R sour cc <C r C a US tic)- 

In this paper, we estimated the microlensing time scales 
for the X-ray, UV and optical emitting regions of the accre- 
tion disc using the following three methods: 

(i) by converting the distance scales of microlensing 
events to the corresponding time scales according to the for- 
mula (13) in which i? SO urcc is replaced by the distance from 
the center of accretion disc. Caustic rise times {ta me) are 
then derived from the simulated variations of the normal- 



Figure 4. The same as in Fig. 2. in the case when time scale 
corresponds to z t = 0.04, z s = 1.69 (i.e. to Q2237+0305 lens 
system) . 

ized total flux in the X-ray, UV and optical spectral bands 
(see Figs. 2-4) by measuring the time from the beginning to 
the peak of the magnification event 

(ii) by calculating the caustic times (t C austic) according to 
equation (14) 

(iii) using light curves (see Figs. 5 and 6) produced when 
the source crosses over a magnification pattern. Rise times 
of high magnification events {tn me) are then measured as 
the time intervals between the beginning and the maximum 
of the corresponding microlensing events (for more details, 
see the next section). 



4 RESULTS AND DISCUSSION 

In order to explore different cases of microlensing and eval- 
uate time scales for different spectral bands, first we numer- 
ically simulate the crossing of a straight-fold caustic with 
parameters Ao = 1, /3 = 1, k = +1 and r caU stic = 9000 R g 
over an accretion disc with an inclination angle 35°, that is 
stratified into three parts: 

(i) The innermost part that emits X-ray continuum 
(1.24 A - 12.4 A or 1-10 keV). The inner radius is taken to 
be Ri„„ = Rms (where R m3 is the radius of the marginally 
stable orbit: R„ l3 — 6 R g in the Schwarzschild metric) and 
outer radius is R ou t = 80 R g (where R g = GM/c 2 is the 
gravitational radius for a black hole with mass M). 

(ii) An UV emitting part of the disc (contribute to the 
emission from 1000 A - 3500A), with R lnn = 100 R g and 
Rout = 1000 R g . 

(iii) An outer optical part of the disc with Ri nn = 100 
R g and R ou t = 2000 R g that emits in the wavelength band 
from 3500 A until 7000 A. 

Having in mind that the aims of this investigation are 
to study the microlensing time scales for different emitting 
regions and time dependent response of amplification in dif- 
ferent spectral bands, we considered microlensing magnifi- 
cation patterns only for image A of Q2237+0305 and for 
a "typical" lens system. Our intention was not to create a 
complete microlensing model for a specific lens system, and 
therefore we did not analyze the differences between images 
(as for instance the time delay between them). The vari- 
ations in the total flux in the different spectral bands are 
given in Figs. 2-6. In Figs. 2 and 3 the simulations for a 
typical lens system with zi — 0.5 and z s — 2 and for two 
different "caustic sizes" are given. As one can see from Figs. 
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Table 1. The estimated time sc ales (in years) fo r microlensing of the X-ray, UV and optical emission region for lensed QSOs observed 
by Chandra X-ray Observatory jDai et al. Il2004) . The calculated caustic times t caU stic are obtained according to formula (14) for the 
following values of the cosmological constants: Ho = 75 km s — iJVlpc -1 and Co = 0.3. The caustic rise times tuME are derived from 
caustic crossing simulations (see Figs 2 - 4). The black hole mass is assumed to be 10 8 Mq. 



Object 






X- 


-ray 


UV 




optical 








^caustic 


tHME 


^-caustic 


tHME 


^caustic 


tHME 


HS 0818+1227 


3.115 


0.39 


0.572 


0.660 


7.147 


7.070 


14.293 


15.160 


RXJ 0911.4+0551 


2.800 


0.77 


0.976 


1.120 


12.200 


12.080 


24.399 


25.880 


LBQS 1009-0252 


2.740 


0.88 


1.077 


1.240 


13.468 


13.330 


26.935 


28.570 


HE 1104-1805 


2.303 


0.73 


0.918 


1.050 


11.479 


11.370 


22.957 


24.350 


PG 1115+080 


1.720 


0.31 


0.451 


0.520 


5.634 


5.570 


11.269 


11.950 


HE 2149-2745 


2.033 


0.50 


0.675 


0.780 


8.436 


8.350 


16.871 


17.890 


Q 2237+0305 


1.695 


0.04 


0.066 


0.080 


0.828 


0.820 


1.655 


1.760 



Table 2. Average rise times (< tnME >) of high magnification events, their average number (< iV caus ti c >o) per unit length (Re) and 
their average number (< -/V caust ; c > y ) per year in the light curves of Q2237+0305A (Fig. 5) and "typical" lens system (Fig. 6). 









Q2237+0305 






"Typical" lens 




Disc path 


Sp. band 


< tnME > 


< -^caustic >0 


< -^caustic >y 


< tHME > 


< ^caustic >0 


< -^caustic >y 




X 


0.37 


0.95 


0.20 


3.94 


1.63 


0.06 


y = o 


UV 


1.29 


0.59 


0.12 


10.96 


0.78 


0.03 




Optical 


2.96 


0.51 


0.10 


15.81 


0.62 


0.02 




X 


0.57 


0.52 


0.11 


2.77 


1.65 


0.06 


y = -x 


UV 


1.92 


0.36 


0.07 


6.62 


0.49 


0.02 




Optical 


3.98 


0.26 


0.05 


17.04 


0.38 


0.01 




X 


0.68 


1.61 


0.33 


4.52 


1.40 


0.05 


x = 


UV 


1.38 


0.88 


0.18 


13.25 


0.54 


0.02 




Optical 


2.96 


0.44 


0.09 


31.62 


0.31 


0.01 



2 and 3, the microlensing time scales are different for dif- 
ferent regions, and the durations of variations in the X-ray 
are on order of several months to a few years, but in the 
UV/optical emission regions they are on order of several 
years. Also, as one can see in Figs. 2 and 3, the time scales 
do not depend on "caustic size" which, on the other hand, 
affects only the maximal amplifications in all three spectral 
bands. The results corresponding to the lens system of QSO 
2237+0305 (z t = 0.04, z s = 1.69) are given in Figs. 4 and 
5. We considered the straight-fold caustic (Fig. 4) and a mi- 
crolensing pattern for the image A of QSO 2237+0305 (Fig. 
5). As one can see from Figs. 4 and 5, a higher amplification 
in the X-ray continuum than in the UV/optical is expected. 
In this case, the corresponding time scales are much shorter 
and they are from a few days up to a few months for X-ray 
and a few years for UV/optical spectral bands. The similar 
conclusion arises when we compare the latter results with 
those for a "typical" lens system in the case of microlensing 
pattern (Fig. 6). 

We also estimated time scales for seven le nsed QSOs 
which have been observed in the X-ray band l)Dai et al. I 
2004). For each spectral band two estimates are made: caus- 
tic time (tcaustic) - obtained from equation (14) and caustic 
rise time (tHME) - obtained from caustic crossing simula- 
tions (see Figs. 2-4). In the second case, we measured the 
time from the beginning of the microlensing event until it 
reaches its maximum (i.e. in the direction from the right to 
the left in Figs. 2-4). The duration of the magnification 



event beyond the maximum of the amplification could not 
be accurately determined because of the asymptotic decrease 
of the magnification curve beyond the peak. Estimated time 
scales for different spectral bands are given in Table 1. As 
one can see from Table 1, the microlensing time scales are 
significantly smaller for the X-ray than for the UV/optical 
bands. 

Unamplified and amplified brightness profiles of the X- 
ray emitting region, corresponding to the highest peak in 
Fig. 5 (top right) are presented in Fig. 7. As one can see in 
Fig. 7, the assumed brightness profile of the source is very 
complex due to applied ray tracing method, which allows 
us to obtain an image of the entire disc, not only its one 
dimensional profile. Therefore, we could not use a simple 
source profile for the calculat i on of microlensing time scales 
(as it was done bv lWitt et alll|l995h ). Instead, we estimated 
the frequency of high amplification events (HMEs), i.e. the 
number of such events per unit length, directly from the 
light curves presented in Figs 5 and 6. For models with non- 
zero shear, this frequency depends on the direction of motion 
and for both calculated maps (Q2237+0305A and "typical" 
lens) we counted the number of high magnification events 
along the following paths (see Fig. 1): i) horizontal (y = 0) 
in the direction from ~x to x, ii) diagonal (y = —x) in 
direction from —x to x and iii) vertical (x — 0) in direc- 
tion from y to — y. For each map the lengths of the hori- 
zontal and vertical paths in the source plane are the same 
(13.636 Re for Q2237+0305A and 12.875 R E for "typical" 
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Figure 5. Variations in the X-ray (solid), UV (dashed) and optical (dotted) spectral bands corresponding to the horizontal (top left), 
diagonal (top right) and vertical (bottom) path in the magnification map of the QSO 2237+0305A (Fig. 1 left). 



lens), as well as are the corresponding crossing times (66.22 
years for Q2237+0305A and 337.32 years for "typical" lens). 
For Q2237+0305A the length of the diagonal path in source 
plane is 19.284 Re and the crossing time is 93.62 years, while 
the corresponding length in the "typical" case is 18.208 Re 
and the crossing time is 477.04 years. 

HMEs are asymmetric peaks in the light curves which 
depend not only on microlens parameters but also on the 
sizes of emitting regions in the following sense: the larger 
emitting regions are expected to produce smoothe r light 
curves and more symmetric peaks (|Witt et al.lll995l ). Con- 
sequently, it can be expected that the majority of HMEs 
should be detected in X-ray light curves, less of them in UV 
and the smallest number in optical light curves. Therefore, 
we isolated only clearly asymmetric peaks in all light curves 
and measured their rise times tnME as the intervals between 
the beginning and the maximum of the corresponding mi- 
crolensing events. In the case of Q2237+0305A we found the 
following number of HMEs: i) horizontal path: 13 in X-ray, 8 
in UV and 7 in optical band, ii) diagonal path: 10 in X-ray, 



7 in UV and 5 in optical band and iii) vertical path: 22 in 
X-ray, 12 in UV and 6 in optical band. In case of "typical" 
lens these numbers are: i) horizontal path: 21 in X-ray, 10 
in UV and 8 in optical band, ii) diagonal path: 30 in X-ray, 
9 in UV and 7 in optical band and iii) vertical path: 18 in 
X-ray, 7 in UV and 4 in optical band. 

The average number of caustic crossings per unit length 
(< A^austic >o) and per year (< A?" caus tic > y ) are given 
in Table 2. This table also contains the average rise times 
(< titME >) in all three spectral bands, derived from the 
rise times {tsME) of individual HMEs which are presented 
in Figs. 8 and 9 in the form of histograms. These results, 
as expected, show that the rise times are the shortest and 
the frequency of caustic crossings is the highest in the X- 
ray spectral band in comparison to the other two spectral 
bands. One can also see from Tables 1 and 2 that in the 
case of Q2237+0305A, the average rise times of HMEs for 
all three spectral bands (obtained from microlens magnifica- 
tion pattern simulations) are longer than both, caustic rise 
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Figure 6. The same as in Fig. 5 but for the magnification map of a "typical" lens system (Fig. 1 right). 




Figure 7. Unamplified (left) and amplified (right) brightness profile of the X-ray emitting region, corresponding to the highest peak 
Fig. 5 (top right). The profiles are obtained using the ray tracing method (see, e.g., Popovic et al. 2003a, b and references therein). 
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Figure 8. Rise times (tHME) of high magnification events for all three spectral bands in the light curves of Q2237+0305A (Fig. 5). Top 
left panel corresponds to the horizontal, top right to the diagonal and bottom to the vertical path of accretion disc. 




Figure 9. The same as in Fig 8. but for the light curves of "typical" lens system (Fig. 6) 
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times (obtained from caustic simulations) and caustic times 
(calculated from equation (14)). 

Microlensing can result in flux anomalies in the sense 
that differe nt image flux ratios are ob s erved in different spec - 
tral bands |Popovic fe Chartaj|200l|Popovic et al.ll2006bl ). 
As shown in Figs 2-6. the amplification in the X-ray band 
is larger and lasts shorter than it does in the UV and opti- 
cal bands. Consequently, monitoring of lensed QSOs in the 
X-ray and UV/optical bands can clarify whether the flux 
anomaly is produced by CDM clouds, massive black holes or 
globular clusters (millilensing) or stars in foreground galaxy 
(microlensing). 



5 CONCLUSION 

In this paper we calculated microlensing time scales of dif- 
ferent emitting regions. Using a model of an accretion disc 
(in the center of lensed QSOs) that emits in the X-ray and 
UV/optical spectral bands, we calculated the variations in 
the continuum flux caused by a straight-fold caustic crossing 
an accretion disc. We also simulated crossings of accretion 
discs over microlensing magnification patterns for the case 
of image A of Q2237+0305 and for a "typical" lens system. 
From these simulations we concluded the following: 

(i) one can expect that the X-ray radiation is more am- 
plified than UV/optical radiation due to microlensing which 
can induce the so called 'flux anomaly' of lensed QSOs. 

(ii) the typical microlensing time scales for the X- 
ray band are on order of several months, while for the 
UV/optical they are on order of several years (although the 
time scales obtained from microlensing magnification pat- 
tern simulations are longer in comparison to those obtained 
from caustic simulations). 

(iii) monitoring of the X-ray emission of lensed QSOs 
can reveal the nature of 'flux anomaly' observed in some 
lensed QSOs. 

All results obtained in this work indicate that monitor- 
ing the X-ray emission of lensed QSOs is useful not only 
to discuss the nature of the 'flux anomaly', but also can be 
used for constraining the size of the emitting region. 
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